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Abstract 

We present a renormalizable non supersymmetric Grand Unified 50(10) model 
which, at the price of a large fine tuning, is compatible with all compelling phenomeno- 
logical requirements below the unification scale and thus realizes a minimal extension 
of the SM, unified in S'0( 10) and describing all known physics below Mqut- These 
requirements include coupling unification at a large enough scale to be compatible with 
the bounds on proton decay; a Yukawa sector in agreement with all the data on quark 
and lepton masses and mixings and with leptogenesis as the origin of the baryon asym- 
metry of the Universe; an axion arising from the Higgs sector of the model, suitable to 
solve the strong CP problem and to account for the observed amount of Dark Matter. 
The above constraints imposed by the data are very stringent and single out a par- 
ticular breaking chain with the Pati-Salam group at an intermediate scale Mj ~ 10^^ 
GeV. 
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1 Introduction 



The first LHC runs at 7-8 TeV have led to the discovery of a candidate Higgs boson and 
to the non observation of new particles or exotic phenomena: no signals of new physics have 
been found neither by direct production of new particles nor in electroweak precision tests 
nor in flavour physics. The Standard Model (SM) has so far been confirmed by experiment 
beyond all expectations. This is surprising since the hierarchy problem [1] and, to some 
extent, the elegant WIMP (Weakly Interacting Massive Particle) solution of Dark Matter 
strongly suggested the presence of new physics near the Fermi scale. But as well known 
the hierarchy problem is one of "naturalness": the SM theory is renormalizable, finite, well 
defined and predictive once the dependence on the cut off is absorbed in a redefinition 
of masses and couplings. Thus the theory can indeed work in practice and the hierarchy 
problem only arises at the conceptual level if one looks at the cut off as a parameterization 
of our ignorance on the new physics that will modify the theory at large energy scales. 

The hierarchy problem is not the only naturalness problem in fundamental physics: the 
observed value of the cosmological constant A also poses a tremendous, unsolved naturalness 
problem [2]. Yet the value of A is close to the Weinberg upper bound for galaxy formation [3]. 
According to the anthropic point of view, possibly our Universe is just one of infinitely 
many bubbles (Multiverse) continuously created from the vacuum by quantum fluctuations. 
Different physics takes place in different Universes according to the multitude of string theory 
solutions (~ 10^^^) Perhaps we live in a very unlikely Universe but the only one that 
allows our existence [5], [6]. In the context of the SM one can argue against this view since 
plenty of models have been formulated that easily reduce the fine tuning from 10^"^ to 10^: so 
why make our Universe so terribly unlikely? If to the SM we add, say, supersymmetry, does 
the Universe become less fit for our existence? By comparison the case of the cosmological 
constant is a lot different because the context is not as fully specified as for the SM 

However, as the criterium of naturalness has so far failed, we are lacking at present a 
reliable argument on where precisely the new physics threshold should be located. Because 
of the serious arguments against applying the anthropic philosophy to the SM, many still 
remain confident that some new physics will appear not too far from the weak scale; still, 
given the failure so far to detect new physics, there has been a revival of models that 
ignore the fine tuning problem while trying to accommodate the known facts. For example, 
several fine tuned Supersymmetric (SUSY) extensions of the SM have been studied like split 
SUSY [7] or large scale SUSY [8l|9]. There have also been reappraisals of non SUSY Grand 
Unified Theories (GUT) where again one completely disregards fine tuning: in some recent 
papers [TUHT?j several aspects of the problem have been discussed. Here we would like to 
establish whether a relatively simple non SUSY 5*0 (10) GUT extension of the SM exists 
which is able to reproduce all the really compelling data that demand new physics beyond 
the SM at scales at and below Mqut- We consider here a renormalizable model but in 
principle one could also study the same problem in 5*0 (10) versions where each large Higgs 
representation is replaced by products of smaller field multiplets. In our resulting model 
the SM spectrum is completed by the just discovered light Higgs and no other new physics 
is present in the LHC range. At the GUT scale of Mqut > 10^^ GeV the unifying group 
is ^O(IO), broken down to the Pati-Salam group SU{A)c ® SU{2)l ® SU{2)r {Ac2l2r for 
shorthand), valid at an intermediate scale, typically Mj ~ 10^° — 10^^ GeV. Note that, in 
general, unification in SU (5) would not work because a group of rank larger than 4 is required 
in order to allow for a two step (at least) breaking needed, in the absence of SUSY, to restore 
coupling unification and to avoid a too fast proton decay (an alternative within SU{5) is 
to assume some ad hoc intermediate threshold with a set of new particles that modify the 
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evolution towards unification [S]). The proposed non-SUSY SO{10) model, with a single 
intermediate breaking scale M/ between Mqut and the electroweak scale, is compatible with 
the following requirements: 

1- unification of couplings at a large enough scale Mqut compatible with the existing 
bounds on the proton life-time; 

2- a Yukawa sector in agreement with all data on flavour physics, fermion masses and 
mixings, also including neutrinos, as well as with leptogenesis as the origin of the 
baryon asymmetry of the Universe; 

3- an axion, which arises from the Higgs sector of the model, suitable to solve the strong 
CP problem and account for the observed amount of Dark Matter. 

There is no item in the list concerning the onset of inflation in that we assume that the 
infiaton is a gravitationally coupled singlet (or even a larger sector of particles) without SM 
gauge interactions. It turns out that imposing all these requirements is very constraining, so 
that most of the possible breaking chains of 5*0(10) must be discarded and the PS symmetry 
at the intermediate scale emerges as the optimal solution. We show that all these different 
phenomena can be satisfied in this fully specified, although schematic, GUT model, with a 
single intermediate scale at Mj ~ 10^^ GeV. In fact, within this breaking chain, the see-saw 
and leptogenesis mechanisms can both be made compatible with Mj ~ 10^^ GeV, which 
is consistent with the theoretical lower limit on the lightest heavy right-handed neutrino 
for sufficient leptogenesis [13] given by Mi > 10^ GeV. The same intermediate scale M/ 
is also suitable for the axion to reproduce the correct Dark Matter abundance. Given the 
indicative character of our study we limit ourselves to Leading Order (LO) evolution, a 
relatively economic choice of SO{10) Higgs multiplets, a crude threshold matching and a 
sketchy approach to leptogenesis and to axions. We are aware that such a GUT model is 
terribly fine tuned, because of its explicit hierarchy problem that is also manifest in the 
necessary huge splittings within the Higgs multiplets (a sort of generalized doublet-triplet 
splitting problem). We consider this model as an extreme reference case where a minimum 
of new physics is introduced to realize SO{10) Grand Unification and to accommodate all 
unavoidable experimental requirements. Note that, given the experimental values of rriH, frit 
and as{Tnz)i the vacuum instability occurring in the SM near the intermediate scale Mi [12] 
is maintained in this model because below Mi the coupling evolution is the same as in the 
SM. However the evolution is somewhat distorted above the intermediate scale. 

Based on the earlier results in [T0]4T3|[T7M2T] . we are led to the following breaking chain: 

Som ''^-I^'^" Ac2L2n ^c^lIy 3c ly , (1) 

with the Pati-Salam (PS) group Aq^lI^ being the intermediate gauge symmetry group. For 
the breaking of 5*0(10) to the PS group we adopt a 210// of Higgs, which gives suitable values 
for Mqut and Mj. The representation 54 also leads to the PS intermediate group but it 
leaves the left-right symmetry unbroken and, moreover, the mass scales Mqut and Mj turn 
out to be not appropriate. As discussed in Sect. 3 different intermediate groups like Sc^l^rI 
or Aq^lIr fail to end up with a running compatible with all the constraints (we confirm and 
extend to our case the results of Ref. [21]). The breaking down to the SM is achieved by 
using a 126^ (and a less relevant 45^) and the final step to the 3c ly is obtained by means 
of a lO/f. In order to have a suitable axion Dark Matter candidate we also introduce the 
45// representation, with a mass scale close to M/ and a specified transformation property 
under a Peccei-Quinn symmetry. Both 126// and 10// are also necessary to generate fermion 
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210h (1 



16 
126h 



(1,2,2) © (6,1,1) 
(4,2,1) © (4,1,2) 
(1,1,3) © (1,3,1) © (6,2,2) © (15,1,1) 
(6,1,1) © (10,1,3) © (10,3^) © (15,2,2) 
1,1)©(15,1,3)©(15,1,1)©(15,3,1)© (10,2,2) ©(10,2,2)©(6,2,2) 



Table 1: Higgs multiplets under the Pati-Salam group. 



masses (for simplicity, we do not introduce a 120/^ that, in principle, could also contribute to 
fermion masses). We perform a full 3-generation study of fermion masses and mixings. The 
previous works that are closest to the present approach are those in refs. p ^ [T2 | IT3 | [2 m [22]. 
We differ from Ref. [10] in the symmetry breaking chain (those authors choose a 54 for the 
breaking at Mqut) and for the additional representation 45 h involved in the PQ symmetry. 
Also they do not make a full 3-generation fit and do not discuss leptogenesis. In Ref. [T2] 
there is a detailed fit of fermion masses but no discussion of leptogenesis, axions and the 
5*0(10) breaking chain. We differ from Ref. [13] in that the authors assume a particular form 
of lepton-quark symmetry, do not discuss the PQ symmetry and do not perform a detailed 
fit of all fermion masses. Finally, in Ref. [20] and Ref. [22] the same Yukawa sector as ours 
is considered but neither the implications for leptogenesis nor for axions are discussed. 

The Dark Matter problem is one of the strongest evidences for new physics. In this model 
it should be solved by axions [231426] . It must be said that axions have the problem that 
their mass should be adjusted to reproduce the observed amount of Dark Matter. In this 
respect the WIMP solution, like the neutralinos in SUSY models, is more attractive because 
the observed amount of Dark Matter is more guaranteed in this case. Neutrino masses 
and mixing originate in this GUT extended SM from lepton number violation, Majorana 
masses and the see-saw mechanism while baryogenesis occurs through leptogenesis [27] . One 
should one day observe proton decay and neutrino-less beta decay. None of the alleged 
indications for new physics at colliders should survive (in particular even the claimed muon 
(g-2) discrepancy [221122] should be attributed, if not to an experimental problem, to an 
underestimate of the theoretical errors or, otherwise, to some specific addition to the model 
[3U]). This model is in line with the non observation of — )■ 67 at MEG [21], of the electric 
dipole moment of the neutron [32] etc. It is a very important challenge to experiment to 
falsify this scenario by establishing a firm evidence of new physics at the LHC or at another 
"low energy" experiment. 

The plan of this article is as follows: in Sectj2]we elucidate the role of the various Higgs 
representations useful for our purposes; in Sectj3] we discuss the evolution of the gauge 
couplings from Mz up to the GUT scale, with the PS intermediate gauge group, and argue 
on why other possible breaking patterns are less suitable to accommodate the requirements 
listed in the above-mentioned points 1-3. Sect. 4 is devoted to the study of the fit of the 
fermion masses, mixing and leptogenesis whereas in Sect. 5 we present the implication of our 
results on the mass and Dark Matter contribution of cold axions. In Sect. 6 we draw our 
conclusions. 

2 Higgs representations 

In the following, we describe in more detail the 5*0(10) representations needed to realize 
the program outlined above. It is useful to classify the various submultiplets in terms of 
their PS quantum numbers; they are reported in TablH 
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In the Yukawa sector, the (1,2,2) of the IOh representation can be decomposed into 
(1,2,2) = (1,2,+|) © (1,2,-i) = Hu®Hd under the 3c2l1y group; if 10^ = 10^^ then 
H*^ = Hd as in the SM. It has been shown, for instance in [10], that in the limit of Vcb = 
(with Vij the CKM matrix) the ratio mt/rrib should be close to 1, thus contradicting the 
experimental fact that, even at the GUT scale, mt/rrib ^ 1. On the other hand, although 
the IOh is a real representation from the SO{10) point of view, its components can be 
chosen either real or complex. In the latter case, 10^^ ^ 10*fj and then H* ^ Hd- An 
extra symmetry, in our case the Peccei-Quinn U{\)pq, is present in the lagrangian to forbid 
the Yukawa couplings related to 10|^ (see below for a detailed discussion). This solves the 
problem and helps in keeping the parameter space at an acceptable level. In the following 
for vacuum expectation values (vevs) we will use a short-hand notation like: 

^M.d = ((1) 2, 2)„^rf)io . (2) 

These vev's, of the order of the EW scale, generate fermion masses and break Sq^lIy down 
to 3c1q. The 126// contains two important vevs: 

w/j= ((10,l,3))i26, = ((15,2,2),,d)i26. (3) 

The first one, vr, of order M/, is needed to generate the right-handed neutrino mass matrix 
and to break the PS group down to the SM; while Vu,d contribute to the fermion mass 
matrices and thus must be of the order of the EW scale. We assume that only a single light 
Higgs doublet remains with components in different PS representations. For this to work, an 
0{1) mixing between the (1,2,2) of the 10// and the (15,2,2) of the 126// is needed in the 
effective electroweak doublet; such a mixing occurs below the PS mass scale, so the (15, 2, 2) 
fields must have masses around Mj. 

The role of the 45// representation is better understood in connection with the U(1)pq 
Peccei-Quinn symmetry [23l|24] and the axion solution of the Dark Matter problem [251126] . 
With a complex 10// we can transform the fields participating to the fermion mass matrices 
as follows: 

16// ^ e*"16// , IOh ^ e^^*"10// , 126// -> e~2^"T26// . (4) 

This solves the problem of the Yukawa couplings of the 10^ discussed above and leads to 
the axion as Dark Matter candidate. No other multiplets in the model have a non vanishing 
PQ charge, except for the 45// (see below). It is expected that the U(1)pq be broken by a 
nonzero (126//) at the scale of SU(2) /^ breaking, otherwise 10// would drive the U(l) breaking 
to give MpQ ^ M^y, which is ruled out by experiments. Note that the (10, 1, 3) component 
of 126// responsible for the PS breaking contains a component (1, 1, 3, —1) under the group 
SU{3)c X SU{2)l X SU{2)r x U{1)b-l- However a single 126// is not enough since a linear 
combination of the U(1)pq, T^r and B — L remains unbroken. This problem is solved by 
introducing another representation that can be either a 16 [18] or a 45// [19] or another 126, 
as suggested in [TD]. We have studied the evolution of couplings in the three cases and found 
that the 45// representation leads to the most suitable value of M/. The 45// contains the 
multiplet (1, 1, 3) under PS with vanishing B — L number. Assigning to such a multiplet a 
PQ charge a' ^ a and a vev at the same scale M/ we can break the previous degeneracy 
and have a viable axion candidate. 

It is clear that all components of the 126//, 10// and 45// not involved in the breaking 
chain and in the Yukawa sector must live at a much higher scale (the GUT scale) in order 
not to be in conflict with the assumed breaking pattern and with the bounds on proton 
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decay. For this large separation of scales we have no better motivation than to invoke the 
extended survival hypothesis (ESH) p!7 t [T8 | [33] which is the assumption that at any scale, the 
only scalar multiplets present are those that develop VEVs at smaller scales. According to 
this assumption, we summarize in Tabj2]the mass scales of the Higgs components. 





21Qh 


12Qh 


45j^ 




Mgut 


all components 


(6, 1,1), (10, 3,1) 


(1,3,1), (6, 2, 2), (15, 1,1) 


(6,1,1) 


Mi 




(10, 1,3), (15, 2, 2) 


(1,1,3) 




EW 








(1,2,2) 



Table 2: Mass scales of the Higgs components. 



3 Coupling evolution 

In this section we show that the Higgs representations discussed in the previous section 
are enough to guarantee a sufficiently large GUT scale Mqut and an intermediate scale 
Mi compatible with the see-saw mechanism for neutrino masses, an acceptable amount of 
leptogenesis and a viable axion as a Dark Matter candidate. As already specified in the 
Introduction, we restrict ourselves to 1-loop accuracy; the evolution equations between two 
generic scales Mi 2 are then given by the standard formulae: 

«r'(M2) = «-^(Mi)-^log^, (5) 

where the coefficients Oj depend on the representations of fermions and scalars lighter than 
M2. In the remaining of this section, we adopt the following short-hand notation for the 
fields relevant for the coupling evolution: 

(1,2,2) = $, (10,1,3) = Ar, (15,2,2) = E, (1,1,3) = (t. (6) 

Among the Higgs fields only $ is involved in the evolution of the 2i and ly couplings between 
the EW mass and the intermediate mass scale M/; from M/ to Mqut we have to take into 
account the contributions of S, and ex: the fields $ and A/j contribute to a/^Ci ^ and 
E to the SU{2)l group (with coupling and /3-function coefficient 03^) and all of them 
to a2R. 

We have six evolution equations (three below and three above Mi) and six unknown to 
fix: the three couplings at Mi, the unified coupling a^McuT) and the scales Mi and Mqut- 
The following input values and matching conditions are imposed 



asciMz) = 0.1176 ±0.002 a2L(M2) = 0.033812 ± 0.000021 

aiy{Mz) = 0.016946 ±0.000006 (7) 

asciMi) = a^c{Mi) a2L{Mi) = a'^^^Mi) a^^{Mi) = ^a-^{Mi) + ^aJ(i(M,) 
a4.c{MGUT) = 012l{,Mgut) a4c{MGUT) = oi'2l{.Mgut) ■ 

The /3- function coefficients for the various gauge groups (computed following [36] and [37] ) 
are listed in Tab|31 

The numerical RGE solutions are shown in FigUl we obtain: 

= (1.3 ± 0.2) ■ 10^^ Gel/ Mgc/t = (1-9 ± 0.6) • lO^'' Gel^ (8) 
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^2b 


-7 


ly 


41 


7 


2 


28 


6 


10 


3 


3 



Table 3: /3-function coefficients. 
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Figure 1: Evolution of the coupling constants from Mz to Mqut- 




with ~ 36.4 (or au ~ 0.027). The stated errors only take into account the propagated 
uncertainties from the SM coupling constants and the Z boson mass. Two-loop effects in 
the matching conditions and possible threshold effects are not included. The jump of ai^ 
at the intermediate scale Mj is due to the fact that the generator of ly is a sum of terms 
from the two different factors SU{2)fi and SU{4)c in the PS gauge group, with matching 
condition as given in eq.©. 

From the order of magnitude estimate of the proton lifetime, given by: 

afjnil 

and by using the values of Mqut and au obtained from our coupling evolution, we get: 

^model _ 5. ^Q36y^ (10) 

to be compared to the present experimental limit [38] : 

r'="P = r(p^e+7r°) > l.Ol-lO^^y. (11) 

We see that the model is confortably consistent with the experimental limit. 

We have studied alternative breaking chains to verify whether they could be compatible 
with the requirements outlined in the Introduction. If instead of using a 210/;- to break 
S'O(IO) down to the PS groups we had adopted a 54/^, the intermediate gauge group would 
have been PS x P, where P is a parity symmetry that enforces the left and right couplings 
to be equal from Mj to Mqut- In this case, the Higgs boson sector must be modified to 
become more left-right symmetric and the values of the mass scales and couplings would 
be changed. Such a breaking chain has been already discussed in [21]: from their analysis 
and the present bounds on the proton life-time this chain is ruled out. We have repeated 
the analysis, adding the (15, 2, 2) multiplet of the 126/^ representation and a 45/^ (or IQh) 
at the intermediate scale (all not included in [2T]) to get a viable axion candidate. We 
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confirm that, with present data on the proton decay, the PS x P predicts a too short 
hfe-time, r ~ 10'^^ y. Other interesting possibihties for the intermediate gauge groups are 
those involving Sq^l^rIx, where X = [B — L)/2. In this case, adding the contributions 
of (1, 2, 2, 0) of 126h and that of (1, 1, 3, 0) of the A5h (or that of the (1, 1, 2, -1/2) of the 
IQh), we always get a too small intermediate scale, Mj ~ 10^ GeV, incompatible with the 
see-saw mechanism and difficult to reconcile with the axion explanation of Dark Matter. If 
the symmetry group also includes the P parity, we get Mj ~ (0.4 — 1) 10^^ GeV (depending 
on whether we use the 16h or 45^^, respectively) and a value of the proton life-time roughly 
one-to-two orders of magnitude smaller than its present experimental bound. We conclude 
that the 5*0(10) breaking chains 3(72i2/jlx, with or without P-parity, are not suitable to 
fulfill all the phenomenological constraints considered in this paper. 

4 Fermion masses, mixings and leptogenesis 

Fermion masses arise from the following Yukawa interactions (for economy of parameters 
we are omitting the possible contribution of a 120 h) : 



Ly = 16{hlOH + f 126 h)IQ. (12) 

where the coupling matrices h and / are complex symmetric matrices. Recalling the vev 
definitions in eqs. [21 [3] the fermion mass matrices of the model have the following form: 

Mu = hku + fvu, Md = hkd + f Vd 

(13) 

Mi;' = hK-3fv^, Mi = hkd-3fvd, Mlf = fvR. 

We can rewrite the previous expressions in a more compact form, suitable for a fit to masses 
and mixing angles: 

M^ = r,{H + sF), Md = H + F 

(14) 

=r^{H- 3sF) Mi = H - 3F M^" = r^^F , 



where H = h kd, F = f Vd, = ku/ kd, s = Vu/r^ Vd and tr = Vd/vR. Since the charged lepton 
masses are known with high accuracy, we prefer to express M^ and the neutrino matrices in 
terms of Md and Mi (the latter with eigenvalues at Mqut given by [39] me=0. 46965 MeV, 
m^=99.1466 MeV and m^=l. 68558 GeV) as follows [I2]: 

f3 + s I- s 

Mu = l^^Md + ^^Mi 

^ r.(^M, + l±^M.), (15) 

Ml^ = !|i(M,-MO. 

With respect to the existing literature, in our fit we also add the additional requirement 
of a quantitatively successful leptogenesis |10] : 

ryg^^^ = (5.7 ± 0.6) X 10~^° (90% CL - deuterium only) . (16) 
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The obvious way to compute the baryon-to-photon number ratio would be to implement the 
Boltzmann equations in the fit and then evaluate tjb following, for example, the prescription 
given in Refs. |1I1I12] . However, this procedure is too long and complicated for the precision 
we are aiming to in this work. Instead, we can use approximate expressions, depending on 
the heavy neutrino spectrum of the theory. For example, since leptogenesis is expected to 
occur at a temperature of the order of M^,^, if 10^ < M,^^ < 10^^ GeV only the r Yukawa 
coupling is in equilibrium and the muon and electron asymmetries are indistinguishable. In 
this case one can adopt a two-flavor approach with only A^i contributing to the asymmetry 
if the right-handed spectrum satisfies the condition (M2,3 — Mi) /Mi ^ 1 or adding the N2 
contribution if (M2 — Mi)/Mi ~ Since the heavy spectrum and the Dirac mass matrix 

are not known a priori, that is before making the fit to the observables at the GUT scale, 
we adopt the following general algorithm: 

1- we assume to work with a given number of flavours and active right-handed neutrinos; 

2- in the fit we implement simplified solutions of the Boltzmann equations (see, for in- 
stance, prj): 

3- after the fit, we check a posteriori that the adopted assumptions in step (1) are correct; 

4- in the case of a positive answer, we use the heavy spectrum and the Dirac mass matrix 
obtained from the fit to solve numerically the Boltzmann equations and get a more 
precise determination oi rjB- 

Obviously, if the assumptions in (1) are not correct, we have to modify our approximate 
formulae and run the chain again. It has been emphasized in p^H3] that in S'O(IO) with 
see-saw and renormalizable Yukawa couplings, there is a strict quark-lepton relation not 
very suitable for implementing the mechanism of baryogenesis through leptogenesis. In fact, 
the hierarchy of the eigenvalues of the Dirac mass matrix M^, implied by the quark-lepton 
connection, implies a strong hierarchy among the heavy right-handed neutrinos via the see- 
saw formula and this, in general, produces a too small lightest right-handed mass to generate 
a successful leptogenesis. For a possible way out we assume the easiest possibility, that is 
that A'^i and N2 are sufficiently close in mass to contribute both to 775 ; we further assume 
to work in the two-flavour regime, since we expect the lightest right-handed neutrinos with 
mass of the order of Mj ~ 10^^ GeV, then in a range where only the r Yukawa coupling is 
in equilibrium. 

To estimate the number of independent parameters we proceed as follows. Eq.f llSp con- 
tains 24 real parameters in h and / and 5 (in principle) complex Yukawas, for a total of 34 
parameters. Working in the basis where the charged leptons are diagonal allows to remove 3 
angles and 6 phases contained in the unitary matrix W such that M^*"^ = W M^ W'^, so we 
are left with 25 parameters. From eq. ffT^ we see that the 7 quantities h. /, ku, kd, Vu, v^, vr 
only appear in the 3 combinations r„, s and r/j, so we have 6 and not 10 parameters, which 
means 21 independent parameters. Moreover, since and appear as overall factors in 
eq.f llSp . they can be taken as real, so the total number of parameters is 19, of which 7 are 
phases: 12 in the down mass matrix, 3 in the charged lepton mass matrix, 2 contained in 
the complex s parameter and one each in and r^j. We do not fit neither the lepton masses 
nor tr so, in total, we have to estimate 15 real parameters (12 in the down mass matrix, 2 
contained in the complex s parameter and the real r„) to fit the 15 observables summarized 
in TabJH and in eq. f|T6|) . The relations in eq.f lTSj) hold at the GUT scale, so we have to 
take into account the Standard Model running of the quark masses and mixing and charged 
leptons masses from low energies up to Mgut'-, for a precise calculation, the contribution of 
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the Higgs states with masses around Mj (see TabJl]) cannot be ignored for the running from 
Mi to Mgut'i however, their effect amounts to a minor correction because the corresponding 
beta function coefficients are not large and the running involves a logarithmic dependence on 
the ratio Mgut/Mj which is quite smaller than Mj/Mz- The values of the fermion masses at 
Mgut reported in TabJUare taken from [39]; for the neutrino parameters {6\j ad r) we used 
the recent results in and for the CKM parameters the ones quoted in |15]. Although 
the last two sets of observables are taken at the lowest mass scale in the running, we do not 
expect in this case sizable corrections in the evolution, as it can be appreciated comparing 
our Tab JH with Table VI in [12]. Notice that for the data with asymmetric la error, we take 
the value in the center of the interval as the best fit and consider a symmetric la interval. 
This simplifying choice has a negligible impact on the results of our fit. 



m„ (MeV) 


0.495 ±0.185 


\Vus\ 


0.2254 ±0.0006 


rud (MeV) 


1.155 ±0.495 


\Vcb\ 


0.04194 ±0.0006 


(MeV) 


22.0 ±7.0 


\Vub\ 


0.00369 ±0.00013 


me (GeV) 


0.235 ±0.035 


J 


(3.16 ±0.1) X 10-5 


rrik (GeV) 


1.00 ±0.04 


sin^ 


0.308 ±0.017 


mt (GeV) 


74.15 ±3.85 


sin^ 


0.3875 ± 0.0225 


r 


0.031 ±0.001 


sin^ ^^3 


0.0241 ±0.0025 



Table 4: Input values at the scale Mqut = 2 x 10^^ GeV. 

To check whether the model is able to reproduce the experimental values of fermion masses, 
mixing and 77^, we perform a analysis using: 

where the Pi denote the theoretical values of the observables and the Oj are the experimental 
values extrapolated to Mqut (with o", being the corresponding la uncertainties). We get a 
reasonably good minimum of the x^, namely, for 15 data points: 

xLn, = 17.4; (18) 
the best fit solutions (and the related pulls) of the observables are given in Tabl3 



ohs. 


fit 


pull 


ohs. 


fit 


pull 


m„(MeV) 


0.49 


0.03 


\Vus\ 


0.225 


0.038 


md(MeV) 


0.78 


0.75 


\Vcb\ 


0.042 


-0.208 


m,(MeV) 


32.5 


-1.50 


\Vub\ 


0.0038 


-0.659 


me(GeV) 


0.287 


-1.49 


J 


3.1 X 10-5 


0.589 


mf,(GeV) 


1.11 


-2.77 


sin^ 


0.318 


0.611 


mt(GeV) 


71.4 


0.70 


sin^ 6^23 


0.353 


-1.548 


r 


0.031 


0.10 


sin^ ^^3 


0.0222 


-0.758 




5.699 X 10-1° 


-0.001 









Table 5: Best fit solutions for the fermion observables at the scale Mqut = 2 ■ 10"*^^ GeV. 

We observe that all the experimental data are reproduced within 3(T, the largest contributions 
coming from the h quark mass and to a lesser extent also from the c and s masses and the 
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atmospheric neutrino mixing angle 623. One could argue that an extra source of error on 
the quark masses can well arise from the distortions on the mass evolution induced by the 
deformed evolution with respect to the SM that occur above Mj. On the other hand, the 
tendency of the atmospheric angle to drift toward small values (compared to the input in 
TabJl]), is mainly due to the stringent requirement of a successful leptogenesi^. In fact, as a 
check, we have redone the fit only including the 14 observables in TabJH We have obtained a 
very good minimum, xLin/dof ~ 0.95, very close to the results shown in [12]. Equipped 
with the neutrino mass matrices as obtained from the best fit parameter values (which 
indeed made the heavy neutrino masses very hierarchical, M^-^ : Mj^^ '■ ~ 1 : 10 : 100), 
we computed the resulting baryon-to-photon number ratio, obtaining rjs ~ —10^^, which is 
wrong in sign and magnitude. 

Note that our fit has no degrees of freedom (15 observables vs 15 parameters): due to 
the high non linearity of the problem there is no perfect matching but the average squared 
deviation is of order 1. We are also in the position to make a list of predictions 

including the light neutrino masses, the heavy right-handed masses, the values of the three 
CP-violating phases (the Dirac 5 and the two Majorana phases (pi^2i extracted according to 
the convention used in [17] ) and the value of the effective mass in the neutrinoless double 
beta decay rate, mee- They are summarized in the following Tabj6l The light neutrino 
spectrum corresponds to the normal hierarchy case. The rate of neutrinoless double beta 
decay is too small to be detected by experiments planned for the near future. The masses of 
the heavy right handed neutrinos are in the range 10^^-10^^ GeV. The light neutrino masses 
m^^ are reproduced from the see-saw formula ~ y^v^ jM^ with v ~ vtiio^ by Yukawa 
couplings of the order of 10~^-10~^ which come out from fitting the fermion masses; the sum 
of neutrino masses TjUiy. ~ 0.065 eV is compatible with the constraints from Cosmology |48] . 
S^i'i ^ (0.1 - 1) eV at la. 



light V masses (eV) 


heavy v masses (VS^^ GeV) 


phases (°) 


rriee (eV) 


.0046 
.0098 
.0504 


1.00 
1.09 
21.4 


6 = 88.6 
01 = -33.2 
ci>2 = 15.7 


5 X 10"^ 



Table 6: Predicted values of the light neutrino masses, of the heavy right-handed masses, of 
the leptonic CP-violating phases and of rriee- 

Some additional comments on rjB are worthwhile: the right-handed neutrino masses in Tabj6] 
confirm that our decision to work in the two-flavour regime is the most appropriate for the 
study of leptogenesis in this model; in addition, given that (M2 — Mi) /Mi ~ 10%, it was 
correct to take into account the contribution of N2 also (we verified that the contribution 
of N3 is irrelevant). The baryon asymmetry can now be computed solving numerically the 
kinetic equations, obtaining: 

r/B~5xlO-^°, (19) 

in agreement with the experimental value. 

■^Notice that it is a general property of type-I see-saw in renormalizablc 5*0(10) theories to favor small 
atmospheric mixing, at least in the 2-family case |46| . 
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5 Axions as Dark Matter particles 



Here we attempt to give an estimate of the axion mass which is crucial for axion inter- 
pretation of Dark Matter. On the other hand, the axion mechanism gives a solution to the 
strong CP problem without need to impose an additional constraint in the fitting procedure. 
The mass can be computed using 



1 

Z2 



fn f^TT 



l + Z Fa 



(20) 



where z = mu/md and Fa is the axion decay constant. For our numerical estimate, we fix 
= 92 MeV and = 138 MeV and take Fa = Mj = 1.3 x 10^^; with the parameter z in the 
interval 0.35 — 0.60, we get ~ (4.3 — 4.7) x 10~^ eV, compatible with astrophysical bounds 
|19]. In the case where inflation occurs after the PQ phase transition, the cosmological energy 
density of cold axions can be estimated from |19] : 



Uah'^ ^ 0.7 



Fa 



1012 GeV 



a 



(21) 



where h is the present-day Hubble expansion parameter and a is the "initial misalignment 
angle", varying in the interval (— vr, vr). In order to estimate whether our axions are able to 
fill the experimental bound [50] : 



nah"^ = 0.1199 ±0.0027 



(22) 



we use eq. (!2T]) with Fa and a as free parameters and plot the region where the bound (122|) 
is satisfied at 3a level. This is shown in fig. ([2]); since Fa should be of the order of Mj, where 
the PQ symmetry is broken, we consider the range Fa = {1/3 Mi,3Mi), enclosed in the 
horizontal lines. As we can see, values of Fa close to Mj are perfectly viable to make the 



IxlO'^F 



8x10' 



6x10' 



4x10" 



2x10" 








Fa = Ml/3 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

a 



Figure 2: The region in the {a, Fa) parameter space where the energy density of cold axions 
saturates eg. [2^ at 3a level. The horizontal lines correspond to Fa = 1/3 Mj and Fa = 3M/. 



axion of our model the relevant component of the cold Dark Matter, at least for values of 

In the case where inflation occurs after the PQ phase transition (the situation we discussed 
above), the axion is subject to quantum fluctuations, which leaves a distinctive imprint on 
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the cosmic microwave background (CMB) spectrum. When the temperature is comparable to 
^QCDi such fluctuations induce non- vanishing density perturbations and then a contribution 
to the total power spectrum of the CMB. The recent Planck results [51] are in tension with 
the axion contribution to the CMB computed for the values of a and Fa of our model. 
However, there are several ways to avoid this constraint, such as, for example, invoking an 
axion decay constant during inflation larger than the present one [52] , or theories where the 
scalar fields are coupled to gravity in a non-minimal way [53], or even scenarios in which 
the QCD coupling becomes large at an intermediate or high energy scale in the very early 
Universe [5^ . 

6 Conclusions 

The SM has passed all sorts of tests at Colliders and the general expectation of discovering 
new physics at the LHC has been so far frustrated. Yet the SM cannot explain a number of 
phenomena like Dark Matter, Baryogenesis and neutrino Majorana masses not to mention 
Dark Energy, inflation and quantum gravity. Due to these problems the horizon of particle 
physics must necessarily be extended up to Mqut and even to Mpianck- A suitable arena for 
enlarging the SM to very large energy scales is provided by GUT's. In particular SO{10) 
cut's are very attractive with the spectacular success of the 16 dimensional representation 
that reproduces the quantum numbers of all the fermions in one generation, including the 
right handed neutrinos. In this work we have studied a non SUSY 5*0(10) CUT that could 
extend the validity of the SM up to Mqut describing fermion masses and mixings, Majorana 
masses, neutrinos with see-saw, Baryogenesis and Dark Matter explained by axions. The 
possibility of accommodating all compelling phenomena that demand new physics below 
Mqut in a non SUSY 5*0(10) model is highly non trivial. In fact, it singles out a particular 
breaking chain with a Pati-Salam symmetry at an intermediate mass scale M/ ~ 10^^ CeV. 
We have shown that a reasonable fit to the data can be obtained in this framework; of course, 
the price to pay is a very large level of fine tuning. As our goal was to establish an existence 
proof rather than a completely realistic model, we believe that the rough approximations 
used for the running and the matching can be considered as adequate for our purposes. In 
general the idea of an 50(10) CUT is very appealing but all its practical realizations are 
clumsy, more so in the non SUSY case because of the hierarchy problem and the need of 
an intermediate symmetry breaking scale in order to obtain a precise coupling unification 
and a large enough unification scale to be compatible with the existing bounds on proton 
decay. We discussed here a renormalizable model. Alternatively one could allow for non 
renormalizable couplings in order to work with smaller Higgs representations. We see no 
obstacle in principle to produce a successful model also in this case, but we leave it for 
further study. In conclusion we have considered worthwhile to study the possibility of an 
all-comprehensive version of non SUSY 50(10) CUT and we have built and compared in 
detail with experiment an example of such a theory. 
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LHCPHENONET and Invisibles. 



A Best fit parameter values 

Here we list the best fit values of the 15 parameters used in our fit procedure. The 12 
elements in Md dbTQ clS follows: 

/ (-.0034, .0004) (-7.7 X 10-^-. 0098) (-.0112, -.0712) 
Md(GeV)= (-7.7 X 10-^-. 0098) (.0108, .0010) (.2162, .0060) 

\ (-.0112, -.0712) (.2162, .0060) (1.062, -.0584) 

The complex parameter s and the real parameter ry are: 

s = (.37, -.079) r„ = 60.03 . (23) 

The value of rf> can be fixed, for example, from the solar mass difference Am^2 — 7.54 x 10"^ 
eV^; it turns out that r^^ = 1.21 x 10^^ 



13 



References 



[1] E. Gildener, Phys. Rev. D 14, 1667 (1976); L. Maiani, in Proc. of the Summer School 
on Particle Physics, Gif-sur-Yvette, 3-7 Sep 1979, Ed. by M. Davier et al., IN2P3, 
Paris, France (1979); M. Veltman, Acta Phys. Polon. B12 (1981) 437; E. Witten, Nucl. 
Phys. B 188 (1981) 513 and Phys. Lett. B 105 (1981) 267 

[2] J. Frieman, M. Turner and D. Huterer, Ann.Rev.Astron.Astrophys. 46 385 (2008), 
larXiv:0803.0982l 

[3] S. Weinberg, Phys. Rev. Lett. 59, 2607 (1987) 

[4] M. R. Douglas, |ArXiv:hep-th70602266| larXiv:1204.6"626l 



[5] N. Arkani-Hamed, S. Dimopoulos, S. Kach ru, jArX iv:hep-ph/0501082 G. Giudice, 
R. Rattazzi, Nucl.Phys. B757, 19 (2006), ArXiv:hep-ph/0606105| 



[6] G. Giudice, in "Perspectives on LHC Physics", ed. by G. Kane and A. Pierce, World 
Sci..l a^iv:0801. 25621 

[7] N. Arkani-Hamed and S. Dimopoulos, JHEP 0506 (2005) 073 |arXiv:hep-th/ 0405159|; 
G. F. G iudice and A. Romanin o, Nucl. Phys. B 699 (2004) 65, Erratum-ibid. B706, 65 
(2005), || arXiv:hep-ph/0 406088|; N. Arkani-H amed, S. Dimopoulos G. F. Giudice and 
A. Romanino, Nucl. Phys. B 709 (2005) 3 |arXiv:hep-ph/0409232]; M. Binger, Phys. 
Rev. D73 (2006) 095001, [ii^ :hep-ph/0408240|; D. S. M. Alves, E. Izaguirre and J. G. 
Wacker. [arXivill 08. 3390t M. E. Cabrera, J. A. Casas and A. Delgado. [arXiv:1108.3867l 

[8] L. J. Hall and Y. Nomura, JHEP 1003 (2010) 076 |arXiv:0910.2"235] . 

[9] G. Giudice and A. Strumia, Nucl. Phys. B858 (2012) 63. larXiv:1108.6077[ 

[10] B. Bajc, A. Melfo , G. Senjanovic and F. Vissani, Phys. Rev. D 73, 055001 (2006) 
| hep-ph/0510139| . 

[11] S. Bertolini, L. Di Luzio and M. Malinsky, Phys. Rev. D 80, 015013 (2009) 
[arXiv: 0903. 40491 [hep-ph]]; S. Bertohni, L. Di Luzio and M. Malinsky, AIP Conf. 
Proc. 1467, 37 (2012) jarXi v: 1205. 563 7l [hep-ph]] and references therein; L. Di Luzio, 
larXiv::IilO:3210l [hep-ph]. 

[12] A. S. Joshipura and K. M. Patel, Phys. Rev. D 83, 095002 (2011) |arXiv:1102.5"T48l 
[hep-ph]]. 

[13] F. Buccella, D. Falcone, C. S. Fong, E. Nardi and G. Ricciardi, Phys. Rev. D 86, 035012 
(2012) |arXiv:1203.0829l [hep-ph]]. 

[14] G.F. Giudice, R. Rattazzi and A. Strumia, Phys.Lett. B715 (2012) 142. larXiv:1204.5465l 

[15] S. Davidson, E. Nardi and Y. Nir, Phys. Rept. 466, 105 (2008) |arXiv:0802.2"962l [hep- 
ph]]. 

[16] G Degrassi et al, JHEP 1208 (2012) 098, larXiv:1205.6497l 
[17] F. del Aguila and L. E. Ibanez, Nucl. Phys. B 177, 60 (1981). 
[18] R. N. Mohapatra and G. Senjanovic, Z. Phys. C 17, 53 (1983). 
[19] R. Holman, G. Lazarides and Q. Shafi, Phys. Rev. D 27 (1983) 995. 



14 



K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 70, 2845 (1993) |h ep-ph/ 92 09215] . 
N. G. Deshpande, E. Keith and R B. Pal, Phys. Rev. D 46, 2261 (1993). 
L. Lavoura, Phys. Rev. D 48, 5440 (1993) |hep-ph/9306297] . 

R.D. Peccei, H. R. Quinn, Phys. Rev. Lett. 38 (1977) 1440; Phys. Rev. D16 (1977) 1791; 
S. Weinberg, Phys. Rev. Lett. 40 (1978) 223; F. Wilczek, Phys. Rev. Lett. 40 (1978) 279. 

R. D. Peccei, in CP Violation, ed. C. Jarlskog; Adv. Ser. Direct. High Energy Phys. 
(Worldscientific Pub. Co., Singapore, 1989) 503 

J. E. Kim, AIP Conf.Proc. 1200 (2010) 83. larXiv:0909.3908l 

J.E. Kim and G. Carosi, Rev. Mod. Phys. 82:557 r2010). E^v:0807.3125l 

W. BuchmuUer, R.D. Peccei and T. Yanagida, Ann. Rev.Nucl.Part.Sci. 55:31 1,2005; 
|ArXiv:hep^h 0502169 ; C. S. Pong, E. Nardi, and A.Riotto, Adv.High Energy Phys. 
2012 (2012) 158303, larXiv:1301.3"062l W. Buchmuller, larXiv:1212.3"554l 

Muon g-2 Collab., G. W. Bennett et al, Phys. Rev. D73, 072003 (2006); B. L. Roberts, 
Chin. Phys. C34, 741 (2010) 

A. Hoecker and W. Marciano, The Muon Anomalous Magnetic Moment, in Particle 
Data Group, J. Beringer et al, Phys. Rev. D86(2012) 010001 

K. Kannike et al, larXiv:1111.255T] 

J. Adam et al, [MEG Collaboration] , larXiv: 1303J754I 

C. A. Baker et al, Phys. Rev. Lett. 97:131801,2006; [ArXiv:hep-ex/0602020| ^A3 
S. Dimopoulos, H. M. Georgi and , Phys. Lett. B 140, 67 (1984). 

C. Amsler et al. [Particle Data Group Collaboration], Phys. Lett. B 667, 1 (2008). 
G. AltareUi, larXiv:1303.6065l 

I. G. Koh and S. Rajpoot, Phys. Lett. B 135 (1984) 397. 

D. R. T. Jones, Phys. Rev. D 25, 581 (1982). 

K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G G 37 (2010) 075021. 

Z. -z. Xing, H. Zhang and S. Zhou, Phys. Rev. D 77, 113016 (2008) |arXiv: 0712. 11191 
[hep-ph]]. 

F. locco, G. Mangano, G. Miele, O. Pisanti and P. D. Serpico, Phys. Rept. 472, 1 (2009) 
[arXiv:0809.0631 [astro-ph]]; C. S. Fong, M. C. Gonzalez-Garcia and E. Nardi, Int. J. 
Mod. Phys. A 26, 3491 (2011) |arXiv: 1107.5312,, [hep-ph]]. 

S. Blanchet and P. Di Bari, JCAP 0703, 018 (2007) | hep-ph/0607330] . 

A. Abada, S. Davidson, A. Ibarra, F. X. Josse-Michaux, M. Losada and A. Riotto, 



arXiv:hep-ph/0605281 



[43] F. Buccella, D. Falcone and L. Oliver, Phys. Rev. D 83, 093013 (2011) |arXiv:1006.5"698l 
[hep-ph]]. 



15 



[44] G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A. M. Rotunno, Phys. 
Rev. D 86, 013012 (2012) |arXiv:1205.5254l [hep-ph]]. 

utfit.org/ 

B. Bajc, G. Senjanovic and F. Vissani, Phys. Rev. D 70, 093002 (2004) 
[hep-ph/0402140j . 

S. Antu sch, J. K ersten, M. Lindner, M. Ratz and M. A. Schmidt, JHEP 0503, 024 
(2005) |hep-ph/0501272] . 

G. L. Fogh, E. Lisi, A. Marrone, A. Melchiorri, A. Palazzo, A. M. Rotunno, P. Serra 
and J. Silk et a/., Phys. Rev. D 78, 033010 (2008) |arXiv:0805.25T7l [hep-ph]]. 



http: //pdg. Ibl.gov/2009/reviews/rpp2009-rev-axions.pdf 



P. A. R. Ade et al. [Planck Collaboration]. larXiv: 1303.50761 [astro-ph.GO]. 
P. A. R. Ade et al. [Planck Collaboration], |M2CiYil303,5082J [astro-ph.GO]. 
A. D. Linde, Phys. Lett. B 259, 38 (1991). 

S. Folkerts, C. Germani and J. Redondo. larXiv: 1304. 72701 [hep-ph]. 
K. S. Jeong and F. Takahashi. larXiv:1304.8l3T] [hep-ph]. 



16 



